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ABSTRACT 

We study several effects that influence the strength of the proton cyclotron and atomic features in the 
thermal spectra of magnetic neutron stars. We show that it is possible for vacuum polarization to strongly 
suppress the spectral lines when the magnetic field B > 10 14 G. For weaker fields (B < 7 x 10 13 G), the 
surface spectrum is unaffected by vacuum polarization; thus the proton cyclotron absorption line can 
have a large equivalent width. Using an approximate calculation of the synthetic spectra, we show that 
variation of magnetic fields over the neutron star surface leads to broadening of the line. The recently 
detected absorption features from the isolated neutron stars RX J1308.6 + 2127, RX J1605.3 + 3249, and 
RX J0720.4 — 3125 can plausibly be explained by the proton cyclotron resonance, with possible blending 
due to atomic transitions, in the atmosphere of the star. 

Subject headings: line: formation — magnetic fields — stars: atmospheres — stars: magnetic fields — 
stars: neutron — X-rays: stars 
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mosphere emission of NSs as a function of magnetic field 
strength. We show that the absorption features observed 
in the three NSs (RX J1308.6 + 2127, RX J1605.3 + 3249, 
RX J0720.4 — 3125) can possibly be explained by the pro- 
ton cyclotron resonance at (unredshifted) energy 

E Bp = — = 0.63 B u keV, (1) 
rripC 

where B = 10 14 I?i4 G is the magnetic field strength. 
(There may be blending of absorption features from radia- 
tive transitions in neutral H atoms). We also show that, 
owing to the vacuum polarization effect, the EW of the 
line can decrease dramatically as the magnetic field enters 
the magnetar regime (B > 10 14 G). 

In §2, we review and elaborate upon the effect of vac- 
uum polarization on radiative transfer in NS atmospheres. 
In §3, we present our numerical results for the atmosphere 
spectra with different magnetic field strengths; these re- 
sults illustrate the suppression of atmosphere line features 
as the field strength increases. We consider synthetic spec- 
tra from NS atmospheres in §4 and discuss our results in 
§5- 

2. VACUUM POLARIZATION EFFECT ON RADIATIVE 
TRANSFER IN NEUTRON STAR ATMOSPHERES 

Before presenting our numerical results of NS atmo- 
spheric spectra for different magnetic field strengths (§3), 
we first review and elaborate upon the key physics of the 
vacuum polarization effect on radiative transfer in strong 
magnetic fields (see Lai & Ho 2002, Ho & Lai 2003, here- 
after LH02, HL03; see especially §2 of Lai & Ho 2003a, 
hereafter LH03). 

2.1. Mode Conversion: Two Magnetic Field Regimes 

Quantum electrodynamics predicts that in strong mag- 
netic fields the vacuum becomes bircfringent. In a mag- 
netized NS atmosphere, both the plasma and vacuum 

Current address: Kavli Institute for Particle Astrophysics and Cosmology, Stanford University/Stanford Linear Accelerator Center, P.O. Box 
20450, Mail Stop 29, Stanford, CA, 94309 



1. INTRODUCTION 

In the last few years, considerable observational re- 
sources (e.g., Chandra and XMM-Newton telescopes) have 
been devoted to the study of thermal emission from iso- 
lated neutron stars (NSs) and, in particular, to the search 
for spectral features in the radiation. For many NSs, the 
spectra are found to be featureless and often well fit by 
a blackbody (e.g., Marshall & Schulz 2002; Burwitz et 
al. 2003; see Pavlov, Zavlin, & Sanwal 2002 for a review). 

Recently, absorption features have been found in the 
thermal emission of several isolated NSs. For example, 
the spectrum of the young NS IE 1207.4 - 5209 (with sur- 
face temperature T g ~ 2 MK) shows features at 0.7 and 
1.4 keV (Sanwal et al. 2002; Mereghetti et al. 2002; Hailey 
& Mori 2002) and possibly at 2.1 and 2.8 keV (Bignami 
et al. 2003). Two NSs, RX J1308.6 + 2127 (= RBS 1223) 
and RX J1605.3 + 3249, belonging to the class of dim, 
radio-quiet isolated NSs (see Haberl 2003), have been ob- 
served to possess broad absorption features in their spec- 
tra: at ~ 0.2 - 0.3 keV for RX J1308.6 + 2127 (Haberl 
et al. 2003a) and at ~ 0.45 keV for RX J1605.3 + 3249 
(van Kerkwijk et al. 2004). These are in contrast to two 
similar dim NSs, RX J1856.5 - 3754 (Pons et al. 2002; 
Burwitz et al. 2003) and RX J0720.4 - 3125 (Paerels et 
al. 2001), which show featureless, blackbody-like spectra 
(however, see Haberl et al. 2003b for the detection of a line 
at 0.27 keV in RX J0720.4 - 3125). It is particularly strik- 
ing that, although RX J0720.4 - 3125, RX J1605.3 + 3249, 
and RX J1308.6 + 2127 have similar effective tempera- 
tures (~ 1 MK), the equivalent widths (EW) of their 
lines are very different (van Kerkwijk 2003), from weak- 
est in RX J0720.4 - 3125 (EW « 40 eV) to stronger 
in RX J1605.3 + 3249 (EW w 80 eV) to the strongest 
in RX J1308.6 + 2127 (EW w 150 eV). Motivated by 
these observations, we study in this paper the strength 
of the proton cyclotron line (and atomic lines) in the at- 
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polarizations contribute to the dielectric property of the 
medium. A "vacuum resonance" arises when these two 
contributions "compensate" each other (Gnedin, Pavlov, 
& Shibanov 1978; Meszaros & Ventura 1979; Pavlov & 
Shibanov 1979; Ventura, Nagel, & Meszaros 1979). For a 
photon of energy E, the vacuum resonance occurs at the 
density 

p v - 0.964 Y- 1 B\ i E\r' 1 g cm" 3 , (2) 
where Y e is the electron fraction, E\ = E/(l kcV), and 
/ = f(B) is a slowly varying function of B and is of or- 
der unity (LH02; HL03). For p > py (where the plasma 
effect dominates the dielectric tensor) and p < py (where 
vacuum polarization dominates), the photon modes (for 
E <C Ese = 1-16 B14 MeV, the electron cyclotron energy) 
are almost linearly polarized: the extraordinary mode (X- 
mode) has its electric field vector perpendicular to the k- 
B plane, while the ordinary mode (O-modc) is polarized 
along the k-B plane (where k specifies the direction of 
photon propagation and B is the unit vector along the 
magnetic field). Near p = py, however, the normal modes 
become circularly polarized as a result of the "cancella- 
tion" of the plasma and vacuum effects. When a photon 
propagates outward in the NS atmosphere, its polariza- 
tion state will evolve adiabatically if the density variation 
is sufficiently gentle. Thus, a X-mode (O-mode) photon 
will be converted into a O-mode (X-mode) as it traverses 
the vacuum resonance. For this conversion to be effective, 
the adiabatic condition must be satisfied: 

E>E ad = 2.55 (/ tantf feB |l - Ui |) 2/3 ( -J^J keV, 

(3) 

where 9^b is the angle between k and B, u, = 
{Eb p /E) 2 , and H p = \dz/dlnp\ is the density scale 
height (evaluated at p — py) along the ray. For an 
ionized hydrogen atmosphere, H p ~ 2kT/(m p g cos 9) = 
1.65 Te/(gi4 cos 6) cm, where T = 10 6 T$ K is the tempera- 
ture, g — 10 14 gi4 cm s -2 is the gravitational acceleration, 
and 9 is the angle between the ray and the surface nor- 
mal. In general, the mode conversion probability is given 

byPcon = l-e^ /2)(B/ ^ d)3 . 

Because the two photon modes have very different opac- 
ities, the vacuum polarization-induced mode conversion 
can significantly affect radiative transfer in magnetized 
atmospheres. When the vacuum polarization effect is ne- 
glected (nv), the decoupling densities of the O-mode and 
X-mode photons (i.e., the densities of their respective pho- 
tospheres) are approximately given by (see LH02) 

Po ,„v « 0.42 T 6 - 1/4 £ 3/2 S- 1/2 g cm- 3 (4) 

px,nv « 486 T- 1/4 El /2 S-V 2 B 14 g cm- 3 , (5) 

where S = 1 — e~ E l kT . Thus the X-mode photons are pro- 
duced in deeper, hotter layers in the atmosphere than the 
O-mode photons. When vacuum polarization is taken into 
account, the decoupling densities can be altered depend- 
ing on the location of the vacuum resonance py relative 
to po,nv and px,nv (see Fig. 1). For "normal" magnetic 
fields, 

B<B,«6.6x 10 13 Te 1/s Ei 1/4 S-V 4 G, 2 (6) 

2 Since B; depends weakly on T, one may use T c £f/(10 6 K) as an est 
here. 



the vacuum resonance lies outside both photospheres 
{py < po,nv < Px,nv); for the magnetar field regime, 
B > Bi, the vacuum resonance lies between these two pho- 
tospheres (po.nv < Pv < Px,nv)- These two field regimes 
correspond to qualitatively different vacuum polarization 
effects on the radiative transfer. 

We can estimate the effective decoupling (photosphere) 
densities of the two photon modes when mode conver- 
sion at the vacuum resonance is taken into account (see, 
e.g., LH02). Consider first the "normal" field regime 
B < Bi. The free- free opacities of the two modes are 
approximately given by (neglecting angle-dependent fac- 
tors) ko — 9.3piT 6 1 ^ 2 E± 3 S cm 2 g _1 and kx — Kou^ 1 , 
where p = p\ g cm -3 and u c = (EBe/E) 2 >• 1. From 
hydrostatic equilibrium, the column density y at density 
p is y = 0.83 piT 6 g cm" 2 (for g 14 = 2). The col- 
umn depth of the X-mode photosphere is determined by 
fy^Kxdy + fy V Kody = 2/3. Thus the photosphere den- 
sity of the X-mode photons (which convert to the O-mode 
upon traversing the vacuum resonance) is 

1/2 

PX^O - [Px,nv + PV(1 - U c )] 

^Px,„v[l- (5/^)1 V2 • (7) 

Similarly, the photosphere density of the O-mode photons 
(which convert to the X-mode at the resonance) is 

pO^X^PCnvtl+^/Bi) 4 ] 172 - (8) 

Thus, when (B/Bi) 4 <C 1, the photospheres of the two 
modes are given by equations (4) and (5), and there is 
no net change in the total emission spectrum. We note, 
however, that the X-ray polarization signals are dramat- 
ically affected by vacuum polarization in this "normal" 
field regime (see Lai & Ho 2003b). 

For the magnetar field regime, we find that the effective 
photosphere densities are given by 

Po^x^Pv[l + (Si/S) 4 ] 1/2 (9) 

PX^O — PO,nv (10) 

Thus for (B/Bi) 4 ^> 1, the O-mode photosphere is un- 
changed by vacuum polarization, while the X-mode pho- 
tons (which carry the bulk of the flux) emerge from the 
vacuum resonance layers. We expect that the total spec- 
trum will be affected by vacuum polarization in this regime 
(see §3). 

2.2. No Mode Conversion 

The numerical results shown in §3 are based on the 
"full mode conversion" ansatz, i.e., when we treat the 
vacuum polarization effect, mode conversion is assumed 
for all photon energies and directions of propagation. In 
reality, mode conversion is complete only when the adia- 
batic condition is satisfied [see eq. (3)]. A proper treat- 
ment of partial mode conversion in radiative transfer re- 
quires solving the Stokes parameters of the radiation field, 
rather than the intensities of the two photon modes (see 
LH03). The current generation of NS atmosphere models 
cannot properly treat partial mode conversion. We note 

: for Tq. They only differ by a factor of a few in the cases considered 
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Fig. 1. — A schematic diagram illustrating how vacuum polarization-induced mode conversion affects the emergent radiation from a mag- 
netized NS atmosphere. The photosphere is defined by where the optical depth (measured from the surface) is 2/3 and is where the photon 
decouples from the matter. The left side applies to the "normal" field regime [B < 7 X 10 13 G; see eq. (6)], in which the vacuum resonance lies 
outside the photospheres of the two modes. The right side applies to the "superstrong" field regime (B > 7 X 10 13 G), in which the vacuum 
resonance lies between the two photospheres. 



that, in the opposite limit (i.e., the complete non-adiabatic 
limit or when mode conversion at the vacuum resonance 
is turned off), the X-mode photon decoupling depth can 
still be affected by vacuum polarization. This is because 
the X-mode opacity (without mode conversion) exhibits a 
sharp rise near the resonance. As a result, the X-mode 
decoupling density is close to py for sufficiently high pho- 
ton energies (see LH02 and HL03). As discussed in §4 of 
HL03, the proton cyclotron line can be suppressed even 
in the non-adiabatic limit, when the magnetic field is in 
the range a few x 10 14 G < B f- 4 / 5 T 6 1/10 < 2 x 10 15 G. 
(This range of magnetic fields is very approximate since 
they depend on the direction of photon propagation and 
redistribution of photon spectral flux occurs in real atmo- 
spheres.) 

By assuming full mode conversion in §3, we are artifi- 
cially exaggerating the vacuum polarization effect on the 
emergent spectrum; despite this, we expect the trend (i.e., 
the magnetic field dependence of the spectral features) de- 
scribed in §3 to be correct (see HL03 for a comparison 
of numerical results based on the "full mode conversion" 
ansatz and the "no conversion" ansatz). Most importantly, 
near the ion cyclotron energy (ui « 1), the adiabatic con- 
dition is always satisfied [see eq. (3); see especially Fig. 6 of 
LH03]. For example, at B = 7 x 10 13 G (E Bp = 0.44 keV), 
the adiabatic condition is satisfied when E > 0.40 keV 
and E > 0.42 keV for 6 kB = 15° and 30°, respectively; at 
B = 10 14 G {E Bp = 0.63 keV), eq. (3) is satisfied when 
E > 0.54 keV and E > 0.58 keV for 6 kB = 15° and 30°, 
respectively. Thus we expect the mode conversion calcu- 
lation to be more accurate. 



3. ATMOSPHERE SPECTRA: DEPENDENCE ON MAGNETIC 
FIELD STRENGTH 

Figures 2-4 depict the spectra of ionized hydrogen atmo- 
spheres for different field strengths, all with T c g — 10 6 K 
and the magnetic field perpendicular to the stellar sur- 
face. These correspond to emission from a local patch on 
the NS surface and are calculated using atmosphere mod- 
els in radiative equilibrium (see Ho & Lai 2001 and HL03 
for details). When the vacuum polarization effect is ne- 
glected (see the curves labeled "no vp"), all spectra exhibit 
a broad absorption feature due to the proton cyclotron res- 
onance (Ho & Lai 2001; cf. Zane et al. 2001; Lloyd 2003). 
When the vacuum polarization effect is included, the spec- 
tra for different field strengths can be very different. For 
B = A x. 10 13 G, vacuum polarization produces no mea- 
surable difference in the spectrum (see Fig. 2), although 
we note that the vacuum effect increases the temperature 
of the atmosphere layers outside the photospheres of both 
photon modes (see HL03). For B = 10 14 G (see Fig. 4), 
which satisfies (B/Bi) 4 ^> 1, the effective X-mode photo- 
sphere is located at the vacuum resonance density where 
mode conversion takes place [po->x — Pv] see eq. (9)]. 
Since the X-mode photons (which carry the bulk of the 
thermal energy) emerge from shallower layers of the at- 
mosphere (compared to the no vacuum polarization case), 
the high-energy spectral tail is softened, and the spec- 
trum is closer to a blackbody (although the spectrum is 
still harder than blackbody because of non-grey opacities). 
Also, for magnetar field strengths, since photons with en- 
ergies around E Bp (both near the cyclotron line center 
and outside the line) decouple from similar depths (i.e. at 
density py) in the atmosphere, the proton cyclotron line 
is suppressed by the vacuum polarization effect (cf. Zane 
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E (keV) 

Fig. 2. — Spectra of fully ionized hydrogen atmospheres with T cf f = 10 6 K and B = 4 X 10 13 G. The thick solid line refers to the atmosphere 
with the vacuum polarization effect included (vp), the dashed line refers to the atmosphere without vacuum polarization (no vp), and the 
dotted line is for a blackbody with T = 10 6 K. The result for a partially ionized atmosphere including vacuum polarization is also shown 
(thin solid line). 
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Fig. 3.— Same as Fig. 2, except for B = 7 x 10 13 G. 
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Fig. 4. — Same as Fig. 2, except from B = 10 14 G. 
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et al. 2001; Lloyd 2003). We see from Fig. 3 that the 
B = 7 x 10 13 G case is intermediate between the "normal" 
field regime and the magnetar field regime. 

For the temperature and magnetic field strengths con- 
sidered in Figs. 2-4, the fraction of neutral H atoms in 
the atmosphere is not negligible, and atomic bound-bound 
and bound-free transitions may also give rise to absorption 
features. However, partially ionized magnetic atmosphere 
models are not as well developed as fully ionized models 
(see Ho et al. 2003). In Figs. 2 and 4, we also show the 
spectra of partially ionized H models for B = 4 x 10 13 G 
and B = 10 14 G, including the effect of vacuum polariza- 
tion (using the "full mode conversion" ansatz); these are 
based on the latest opacity and equation of state results for 
magnetized H plasmas (Potekhin & Chabrier 2003, 2004). 
For B = 4 x 10 13 G, the absorption feature associated with 
the ground-state bound-free transition lies at ~ 0.44 keV, 
and the bound-bound transition (from s — 0,is — to 
s = 1, v = 0) is at 0.35 keV, which is close to the proton 
cyclotron resonance; these features are unaffected by the 
vacuum polarization effect. For B — 10 14 G, the bound- 
free feature is at 0.54 keV, and the bound-bound transi- 
tion is at 0.74 keV, which is also close to Eb p - Like the 
proton cyclotron resonance in the magnetar field regime, 
these atomic features are also suppressed by the vacuum 
polarization effect. 

The results shown in Figs. 2-4 assume the adiabatic con- 
dition is always satisfied [see eq. (3) and §2.2]. Numerical 
models, with B < 10 14 G and T off = 10 6 K and assum- 
ing no mode conversion, do not show suppression of the 
proton cyclotron line since there is negligible photon flux 
at high energies to be redistributed to lower energies; line 
suppression (with no mode conversion) starts becoming 
important at T cB « 3 x 10 6 K for B « 10 14 G (see, e.g., 
Fig. 8 of Ho et al. 2003). Nevertheless, as discussed in 
§2.2, we believe the full mode conversion calculation to be 



more accurate (since the adiabatic condition is always sat- 
isfied) near the proton cyclotron energy, and the vacuum 
polarization-induced line suppression shown here should 



4. SYNTHETIC SPECTRA FROM NEUTRON STAR SURFACES 

The spectra presented in §3 correspond to emission from 
a local patch of the NS surface. The proton cyclotron line 
in the local spectrum exhibits a rather spiky feature, which 
is clearly not observed (see also Ho & Lai 2001; Zane et 
al. 2001; Lloyd 2003). To compare with observations, we 
need to calculate the synthetic spectrum from the whole 
NS surface. Such synthetic spectra is necessarily model- 
dependent (see Zane et al. 2001 for a different model and 
method), as the variations of magnetic field and tempera- 
ture on the NS surface are unknown. 

Here, for illustrative purposes, we assume a dipole mag- 
netic field configuration and consider two models for the 
temperature variation: (1) a uniform surface tempera- 
ture over the entire NS surface 2b(r); (2) a (somewhat 
arbitrarily-chosen) non-uniform temperature distribution 
given by 

T(r, 7 ) = T p (r) cos 2 7 + T cq (r) sin 2 7> (11) 

where 7 is the magnetic colatitude, r = — J pkq s dr is 
the Thomson depth in the atmosphere, and Kq s is the 
Thomson opacity. We generate the temperature profiles 
T p (r) and T eq (r) from our atmosphere models for effec- 
tive temperatures T c r and magnetic fields B given by 
{T cS ,B) = (T c p ff ,S p ) and {T^,B eq ), respectively, and we 
assume Tq(t) is the same as T p (t). 

Let the angle between the line of sight toward the ob- 
server (the Z-axis) and the axis of the magnetic pole be 
O. (For a spinning NS, obviously varies as the star ro- 
tates.) The observed spectral flux can be calculated using 
the standard procedure, including general relativistic ef- 
fects (Pechenick, Ftaclas, & Cohen 1983; see also Pavlov 
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& Zavlin 2000). At a given local patch on the NS sur- 
face, we denote the direction of the photons which enter 
the detector by k; the angle between k and the local sur- 
face normal is 6 and the azimuthal angle of the emitting 
patch is <j). For a given mode j, we approximate the spe- 
cific intensity of the radiation from this patch by B„(T)/2 
(neglecting the correction to the source function due to 
scattering), where B V (T) is the Planck function and T is 
the temperature at tI{6,4>) = 2/3. Here, ri(9,4>) is the 
effective optical depth along the ray, 

where Kj and are the absorption and total opacities, 
respectively. Interpolating from a table of calculated op- 
tical depths r^(r, 9, </>), we determine the value of r at 
which t 3 v = 2/3, thereby obtaining the source function for 
each mode. The observed flux spectrum /„ (at the ob- 
served frequency vq) from the entire NS surface is then 

/, =/*+/£, with 

fi = (l + z)- 1 (J)' fjdep £ /2 d6 sintfcostf^M, 

(13) 

where (1 + z) = (1 - 2GM/i?c 2 )~ 1 / 2 is the redshift fac- 
tor, 1/q = v/(l + z) and D is the distance to the source 
(we assume the NS mass and radius arc M = 1.4M Q and 
R = 10 km, respectively). 

Figure 5 shows the computed synthetic spectra from a 
fully ionized hydrogen atmosphere with a dipole magnetic 
field distribution (B p = 4 x 10 13 G) and uniform and 
non-uniform temperature distributions over the NS sur- 
face. The proton cyclotron line at Eb p = 0.25 keV for 
B = 4 x 10 13 G is broadened significantly when the emis- 
sion is integrated over the entire NS surface due to the 
variation of the magnetic field: the equivalent widths for 
the two models shown are ~ 80 and 120 eV, as compared 
to ~ 130 and 150 eV for models with a constant magnetic 
field, and the centroid energy shifts by ~ 10—20% (roughly 
similar to the low B results of Zane et al. 2001, despite the 
differences in model and method) . The equivalent width of 
the line is significant, comparable to that of the absorption 
features observed in the isolated NSs RX J1308.6 + 2127 
(Haberl et al. 2003a) and RX J1605. 3+3249 (van Kerkwijk 
et al. 2004) (see §5). 

5. DISCUSSION 

In several recent papers (Lai & Ho 2002, 2003; Ho 
& Lai 2003), we have shown that at superstrong mag- 
netic fields, B > 10 14 G, vacuum polarization can sig- 
nificantly affect the radiation spectrum from magnetized 
neutron star atmospheres: it softens the high-energy tail of 
the spectrum and suppresses the proton cyclotron feature 
(above ~ 0.6 keV, unrcdshiftcd) and other features (see 
also Ho et al. 2003). The latter could provide an explana- 
tion for the non-detection thus far of lines in the observed 
thermal spectra of several anomalous X-ray pulsars (Patel 
et al. 2001, 2003; Juctt et al. 2002; Tiengo et al. 2002; 
Morii et al. 2003) and soft gamma-ray repeaters (Kulkarni 
et al. 2003), which are thought to possess B > 10 14 G. 

3 dc Vries et al. (2004) find long-term changes in the spectra of RX 
of the line. 



In this paper, we have studied the dependence of the 
neutron star atmosphere spectrum on the magnetic field 
strength. As we explain qualitatively and show numeri- 
cally, at normal neutron star magnetic fields, B < 10 14 G, 
vacuum polarization has little effect on the atmosphere 
emission spectra. Therefore, strong proton cyclotron or 
other atomic features may be present in the thermal spec- 
trum. 

Our calculations of neutron star synthetic spectra, tak- 
ing into account the line broadening effect due to magnetic 
field variation over the neutron star surface, show that the 
recently observed broad absorption features in several dim 
isolated neutron stars could be explained naturally as the 
proton cyclotron line (with possible blending from atomic 
lines of neutral hydrogen) from neutron star atmospheres 
with B < 10 14 G (see Fig. 5). The variation in the strength 
of the observed spectral features in these sources is then 
due to different fractions of the surface with B < 10 14 G. 

For RX J0720.4 — 3125, the weakness (phase-averaged 
EW « 40 eV) of the line at 0.27 keV suggests that 
most of the observable surface of this neutron star has 
B > 10 14 G, and the line is from the magnetic equato- 
rial region of the star (where B is weaker) . The observed 
line energy does not change with phase but has a larger 
EW at pulse decline/minimum (Haberl et al. 2003b), 
which also indicates the line is produced at a small re- 
gion near the magnetic equator. Note, however, that Ka- 
plan et al. (2002) and Zane et al. (2002) place an upper 
limit of ~ a few x 10~ 13 s s~ 4 on the period derivative of 
RX J0720.4 - 3125, which, given its period of 8.39 s, im- 
plies a dipole field with -Bdipoio < 5 x 10 13 G; thus our 
inference of the surface field based on the thermal spec- 
trum implies an appreciable non-dipolar magnetic field on 
the neutron star. 

For RX J1308. 6 + 2127, the stronger (EW « 155 eV) line 
at 0.3 keV suggests that most of the observable surface of 
this neutron star has B < 10 14 G. This neutron star also 
has a double peaked pulse profile (Haberl et al. 2003a), 
as compared to the very sinusoidal single peaked pulse 
profile of RX J0720.4 - 3125 (Haberl et al. 2003b). The 
difference in pulse profiles may be due to differences in the 
viewing geometry and angle between the magnetic and 
rotation axes: in the case of RX J1308.6 + 2127, the ob- 
server sees the equatorial region most of the time (thus the 
stronger line and slightly larger pulse fraction) and each 
magnetic pole, while only one pole is visible in the case of 
RX J0720.4 - 3125. 

Our calculations show that, for a given line-emitting 
area, the 0.45 keV line of RX J1605.3 + 3249 should have 
a smaller EW than the 0.27 keV line of RX J0720.4 - 
3125. However, the line is stronger (EW « 80 eV) in 
RX J1605.3 + 3249. One possible way to reconcile the ob- 
served smaller EW of the 0.27 keV line in RX J0720.4 - 
3125 3 is for the line-emitting area (with B < 10 14 G) to be 
a small fraction of the observed surface, while the larger 
fraction of the surface has a much higher magnetic field 
(B > 10 14 G). Alternatively, we have only considered hy- 
drogen, and some of the lines may be due to other ele- 
ments, such as helium. 

As discussed in §2.2, the numerical treatment of the vac- 

)A — 3125; it is not clear whether this variability can affect the EW 
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E (keV) 

Fig. 5. — Synthetic spectra (unredshifted and with R 2 /D 2 factored out) for fully ionized hydrogen atmosphere emission from a NS with 
a dipole magnetic field distribution and uniform and non-uniform [see eq. (11)] temperature distributions over the NS surface. The uniform 
temperature spectrum uses the temperature profile constructed from the atmosphere model with B = 4 X 10 13 G and T c g = 10 6 K, while the 
non-uniform temperature spectrum uses the temperature profiles T p and T eq constructed from the atmosphere models with B = 4 X 10 13 G 
and T c ff = 10 6 K and B = 2 X 10 13 G and T c fj = 5 X 10 5 K, respectively. is the angle between the observer's line of sight and the magnetic 
axis. 



uum polarization effect in the neutron star atmosphere 
models is approximate, i.e., it does not properly account 
for the partial mode conversion associated with the vac- 
uum resonance. Future work in this direction is neces- 
sary (see LH03), as well as a more accurate and com- 
prehensive study of partially ionized atmosphere models 
(see Ho et al. 2003). Furthermore, we have only exam- 
ined thermal emission; there have been cyclotron features 
seen in the non-thermal emission from magnetars (e.g., 
the 5.0 keV feature from SGR 1806 - 20 during out- 
burst, Ibrahim et al. 2002, 2003; the 8.1 keV feature from 
AXP 1RXS J170849 - 400910, Rea et al. 2003). Such 
theoretical/numerical studies, when combined with obser- 
vational data, should provide useful constraints on the na- 



ture of various types of neutron stars, including the dim 
isolated neutron stars. 
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